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Most Markovian queueing networks that arise as models of stochastic 
congestion systems (e.g., communication networks and multiprogrammed 
computer systems) do not have a product form in their stationary probability 
distributions, and hence are not amenable to the simplicity of product-form 
analysis. In this paper we suggest an approach for systematically examining 
the validity of a class of approximation schemes that is based on the idea of 
equivalent networks and is used for the approximate equilibrium analysis of 
nonproduct-form networks. We study equivalent networks, and prove a gen- 
eralization of the so-called "Norton's" Theorem for closed product-form net- 
works in order to study and generalize the equivalent flow method for the 
approximate analysis of nonproduct-form queueing networks. We then present 
the results of a study of the approximation scheme as applied to a type of 
network model (called a central-server model) that arises frequently in mod- 
eling multiprogrammed computer systems. In this model the central server 
uses a priority discipline, so the resulting network is nonproduct form. This 
study demonstrates the situations under which the approximation can be 
expected to do well or poorly and the kinds of errors it introduces. 

I. INTRODUCTION 

Mathematical modeling of stochastic systems frequently gives rise 
to models in a class referred to as Markovian queueing networks — 
specifically, queueing networks whose time evolution can be described 
by a discrete-state, regular Markov stochastic process. Markovian 
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Fig. 1— Notion of an equivalent network; (a) original network, (b) with arrows 
indicating flows between Qi and its complements, and (c) with arrows indicating models 
of flows between Qi and its complement. 



queueing network models, known as product-form networks, have been 
widely studied, owing primarily to their well-understood stochastic 
behavior, and the simplicity of their analysis in equilibrium. However, 
the class of product-form queueing network models is far from ade- 
quate for modeling many simple real-world congestion systems. The 
exact equilibrium analysis of nonproduct-form queueing networks is, 
in most cases, computationally, and often fundamentally, intractable. 
Much effort has, therefore, been directed towards devising approxi- 
mation schemes that attempt to reconcile the conflicting requirements 
of modeling fidelity and the simplicity of product-form analysis. One 
such class of approximation schemes is based on the idea of equivalent 
networks. In this paper we systematically study this approximation. 

By an equivalent network we mean the following (cf. Fig. 1). Con- 
sider a closed queueing network Q constructed from the set of nodes 
M and the subnetwork Qi consisting of nodes M x C M. Let Qi*, be a 
network constructed from Mi such that the joint equilibrium (proba- 
bility) distribution of Qi«, is the same as the marginal joint equilibrium 
distribution of Qi in Q. The network Qieq is then said to be equivalent 
to Qi* Clearly, to study Qi in isolation, one needs to account for the 



* This notion of equivalence may appear unduly restrictive. Why not establish a 
more detailed stochastic equivalence? For the calculation of many performance analysis 
criteria, the present notion is adequate. However, it is easy to see that the equivalent 
networks we later identify yield equality in distribution for the entire process in 
equilibrium. 
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influence of the nodes in M — Mi on the nodes in Q 1# When Q is 
product form, the influence of the complementary network on Q x takes 
an especially simple form, and can be determined by analyzing a 
modified version of the complementary network in isolation! For the 
case where Mi consists of a single node, this fact was first recognized 
by Chandy et al., 1 who called the equivalent network so obtained a 
"Norton" equivalent, because of the similarity of this equivalence to 
Norton equivalence in electrical circuits. 

In Section II we study equivalent networks and demonstrate the 
simplifications that arise for product-form networks. The development 
yields a generalization of Norton's Theorem to multinode subnetworks 
of closed product-form networks. Essentially, the same extension to 
the entire class of closed product-form networks has been obtained 
independently and concurrently by Kritzinger et al. 2 and Balsamo et 
al., 3 through an approach based on verification via detailed computa- 
tions from the product-form solution. Our approach is substantially 
different, in that it derives Norton's Theorem directly as a special case 
of a general result for stochastically equivalent networks. This ap- 
proach is concise, conceptually and intuitively appealing, gives the 
result a probabilistic interpretation, and shows up clearly the role 
played by the product-form solution. It also seems to be the natural 
approach for the purposes of this study. 

This generalization of Norton's Theorem motivates the following 
approximation scheme. Suppose now that Q is a nonproduct-form 
network, but for the purposes of studying the subnetwork Q u we follow 
the equivalence procedure for product-form networks. Suppose also 
that in doing so we find that the version of the complementary network 
we have to analyze, in order to determine the latter's influence on Q u 
is product form. Let the equivalent network thus obtained be $ieq. 
The approximation scheme, referred to above, approximates the equi- 
librium distribution of Qi with that of Qieq (i.e., approximates Q ieq by 
(Jieq). The effort to determine and analyze ietl will, in general, be 
considerably less than the effort to exactly analyze Qi in Q. 

This approximation scheme is an extension of one (often referred 
to in the literature as an equivalent flow approximation) that has been 
utilized by several workers, in the field of network performance anal- 
ysis, with remarkably accurate results. Sauer and Chandy, 4 and Chow 
and Yu 5 use this idea as the basic step in iterative schemes for 
approximating central-server models in which the central server is not 
of product-form type. Schwartz 6 uses the basic scheme directly to 
approximately analyze a model for a multiple-access communication 
system. In Section III we draw upon the theoretical development in 
Section II to study the validity of the approximation scheme when it 
is applied to a simple test-bed model. 
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II. EQUIVALENT NETWORKS 

Consider a closed Markovian queueing network Q consisting of M 
congestion nodes. In this section we study the problem of the equilib- 
rium analysis of a subnetwork Qi (embedded in Q). To simplify the 
discussion we shall limit our considerations to networks of First In, 
First Out (FIFO) nodes. It is easily recognized that the ideas in this 
section can be extended to apply to more general networks. In Section 
2.3 we establish Theorem 1, which explicates the structure of equiva- 
lent subnetworks of the networks described in Section 2.1. By combin- 
ing this result with Theorem 2, we get a generalization of Norton's 
Theorem. 

2. 1 Network specifications 

Q is a closed queueing network consisting of M FIFO nodes (indexed 
by i G {1, • • • , M }). There are R classes/types of customers (indexed 
by r 6 R = {1, •••,R}) with N r customers in the rth class. The M x 
M matrix P (r) = [pj, r) ] is the routing probability matrix of type r 
customers; customers do not change class as they move from node to 
node. For each r in {1, • • • , R], P (r) is a stochastic matrix which, when 
considered as a transition probability matrix, leads to a Markov chain, 
on the state space {1, • • • , M}, with a single positive, communicating 
class. 

Throughout the following discussion, the network state process is 
assumed to be in equilibrium. The state of the ith node (denoted by 
S l ) is a finite string drawn from the set R. Given a state vector S\ r 
G R appearing in the fcth position in the string S' denotes that a 
customer of type r is in the fcth position, in FIFO order, at the node i. 
Thus, by definition, the customer in service is in the first position. S 
= (S 1 , • • ■ , S M ) denotes the state of the entire network Q. The ith 
node is equipped with an exponential server which, when the state of 
the network is S, serves a customer of class r at the rate v ir (S). 

Qi is a subnetwork of Q, consisting of M x (<N) nodes (indexed by 
i G jl, ••• , Mi}). Q 2 is the complementary network consisting of 
M 2 = Id - Mi nodes (indexed by i G {Mi + 1, • • • , M}). 

Some additional notation is inevitable; this we proceed to describe 
in the next subsection. 

2.2 Notation 

N = (iVi, • • • , Nr) is the population vector of the network Q, where 
N r is the number of customers of class r, r G R. 

Let R* = (Un»i{l, • • • , #}") U where denotes the empty string. 
For any s G R*, denote by N r (s) the population of class r in the string 
s, and let 

N(s) = (iVi(s), • • • , N r (s), • • • , N R (s)). 
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For K, a positive integer, let 

Sn = |(*\ • • • , s K ): (for every i, 1 « i < K , s' £ R*) 

and X 7W) = N}. 

As stated in Section 2.1, S = (S\ ■■■ , S M ) denotes the Q-network 
state. Let Si - (S 1 , • • • , S Ml ) denote the Qi-network state and S 2 = 
(S Ml+ \ ■■■ , S M ) denote the Q 2 -network state. Let Fg, and Fg 1 and 
Fg z denote respectively the sets of feasible states, in equilibrium, of 
the state process of the networks Q, Qi , and Q 2 , respectively. 

A network Qi eq constructed from the nodes {1, • • • , Mi} is said to be 
equivalent to Q x if the joint equilibrium (probability) distribution of 
the state processes of Qi eq is the same as the marginal joint equilibrium 
distribution of the state process of Qi in Q. 

2.3 Construction of Qu q 

Let 7r: Fg — > (0, 1) be the equilibrium distribution of the state 
process of the network Q. Let (for every (1 ^ i ^ Af)(l ^ r ^ R)) (for 
every S G Fg) u ir (S) = uUSJ and (for every S, G Fg») 

pf r } A ir \A customer of type r is in service at node 
i/Q\ i is in state Si). 

Construct a network Qi eq from the nodes jl, • • • , Mi} as follows: 

1. The routing between the nodes in Qi eq is the same as in Qi (self 
loops around nodes in Qi are included in Qi eq ). 

2. When the state of Qi eq is S x , node jf(l ^ j ^ Mi) receives an 
exogenous arrival stream of class r customers with (state dependent) 
rate %* Mi+1 pf r ^ lr (S,)p|, r) . 

3. A customer of class r, after completing service at node i (1 < i < 
Mi), leaves the network Qi eq with probability £ JIa^+i P,> r) . 

Theorem 1: Q leq as constructed above is equivalent to the subnetwork 

QiofQ. 

Proof: The intuitive appeal of the construction is manifest. In Step 2 
of the construction, for every i(M x + 1 < i «S M), pp i> ir (Si) is the 
conditional throughput of type r customers through node i, when Q, 
is in the state Si. A fraction pj of this flow through i finds its way 
into node j of Qj. 

A simple detailed proof can be obtained by summing the Kolmogorov 
equilibrium equations for Q over the set \S G Fg : Si fixed}, and 
observing that the resulting equations are exactly the equilibrium 
equations for Qi eq described above (cf. Ref. 7). D 

Remarks: But for the explosion in notation that occurs in setting up 
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a detailed proof, it is clear that the construction of Qi*, described 
above extends easily to networks other than those described in Section 
2.1. In this work, however, we continue to restrict our attention to 
networks of the latter type. 

We now turn to the subclass of product-form networks of the class 
of networks described in Section 2.1. Since we are concerned here with 
FIFO nodes, the service rates cannot be class dependent. We further 
assume that the service rates are not state dependent in any way, i.e., 
we now have 

(for every (1 *£ i «£ M)(l s= r < R) and S G Fg)MS) = "f« 

Let (for every r G R) C (r) C {1, ■ • • , M} be the subset of nodes of Q 
that communicate under P (r) (i.e., in queueing-network terminology, 
the chain corresponding to class r). Let R 2 = jr G R : C (r) n {Mi + 1, 
■ • • , M\ ^ 0) be the set of customer types that visit Q 2 . Let || R 2 II = 
R 2 , || R - R 2 || = Ri (where || || denotes set cardinality), and reindex 
R so that the elements of R 2 receive the highest indices. Let N 2 = 
(Nfl 1+ i, • • • , N R ) and if s is a string in R*, let N 2 (s) = (N Rl+ i(s), ■■■, 
N R (s)), i.e., N 2 (s) is the population vector, of the string s, restricted to 
the classes in R 2 . 

For every N' = (N' Rl+u ■ ■ • , N«) ^ N 2 , consider the network Q 2 (N') 
obtained from Q by replacing all servers in Qi with infinite speed 
servers (i.e., by short-circuiting the nodes in Qi), and placing N' 
customers in the resulting network. Let 7r N - be the equilibrium distri- 
bution of the state process of the network Q 2 (N'). Define for every 
(Mx + Ui$M) ) r6R 2 , 

£*' 4 7T N ' {A customer of type r is in service at node i in Q 2 (N')}. 

2.4 The product-form case 

Theorem 2: If Q is a product-form network then 

for every (Mi + 1 < i « M)(r G R 2 ) and every S t G F& 

pP = £|f-" 2<Sl) (pf r > as defined earlier). 

(Note: it is obvious that (for every (Mi + 1 < i ^ M), r £ R 2 and Si G 
F&)pf;> = 0.) 

Proof: The proof utilizes a simple lemma and is outlined in the 
appendix. D 

Remarks: Theorem 2, when combined with Theorem 1, yields a gen- 
eralization of Norton's Theorem 1 to multinode subnetworks. Even 
though the previous development is specific to the class of networks 
described in Section 2.1, it is clear that the same approach can be used 
to extend Norton's Theorem to the entire class of closed product-form 
networks. The product-form solution continues to play the same role 
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as it does in Theorem 2, i.e., it allows the rates of the external arrival 
streams in Qieq to be computed from an analysis of Q 2 for all possible 
customer populations in Q 2 . 

III. AN APPROXIMATION SCHEME 

In an IBM Research Report, Chow and Yu 6 suggest a somewhat ad 
hoc, iterative approximation scheme for a class of central-server 
models, with a priority discipline at the central server. As mentioned 
earlier, the scheme relies on an inexact application of Norton's Theo- 
rem to such networks. In Section I we described a natural generaliza- 
tion of this so-called equivalent flow approximation scheme to more 
general nonproduct-form networks. In this section, we present the 
results of a detailed study of the application of this approximation to 
a simple, test-bed, central-server network. 

3. 1 The test-bed model 

Consider the two-node network Q shown in Fig. 2. There are two 
customer classes, namely 1 and 2, with Ni and N 2 customers, respec- 
tively (i.e., N = (N U N 2 )). At node 1, the customers of class 1 (high- 
priority) have preemptive priority over class 2 (low-priority) cus- 
tomers; after being preempted by a class 1 customer, when a class 2 
customer reaches the service station again, it resumes service where it 
left off; class 1 and 2 customers have exponential service times with 
rates v\\ and v\ 2 , respectively. Such a service discipline is commonly 
referred to as a preemptive resume discipline. At node 2, there is no 
priority; customers are served in the order in which they arrive 
(irrespective of class), at the class independent exponential service 
rate i> 2 . Customers alternately seek service at nodes 1 and 2 and stay 
in the network forever. This model belongs to a class of central-server 
networks that arise as models of computer systems. 
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3.2 Approximating the test-bed network 

The network described in Section 3.1 is nonproduct form because 
of the preemptive resume discipline at node 1. In order to approximate 
the equilibrium behavior of node 1, we first increase the service rates 
at node 1 to infinity, thus effectively short circuiting the node. Denote 
the resulting network by Q 2 (Fig. 3). Then for each {k u k 2 ) < (Ni,N 2 ) 
analyze Q 2 with ki and k 2 customers of types 1 and 2, respectively, in 
the network. Let (cf. Thm. 2) (for every (k lf k 2 ) < (N U N 2 )) (for every 
r £ {1, 2}) £ 2 k r lM) = Prob {A customer of type r is in service at node 2 
when (k u k 2 ) customers are in Q 2 \. 

This probability will not depend on the sequence in which the (k u k 2 ) 
customers are placed in Q 2 . Since the service rate at node 2 is class 
independent, it is clear that, in equilibrium, all possible states, for any 
arrangement of the customers, are equally likely. From this we can 
directly conclude that 

(for every (0, 0) < (k u k 2 ) < (N lt N 2 )) (for every r £ {1, 2}), 

^ k,+k 2 

Now consider the open network Q Xeq consisting of the node 1 in 
isolation. The service rates and discipline remain the same as in Q. 
When there are n x customers of type 1 and n 2 customers of type 2 in 
(Ji eq then customers of type r (G {1, 2}) enter the network at the rate 
X<"i."2) w here 

(for every (m,rh) *£ (NtM)^^ = ^"""'V 
When a customer finishes service in Q leq , it leaves the network (Fig. 

4)- 
The evolution of the network Q leq can be described by a regular 

Markov process on the state space \(n\,n 2 ) : (n u n 2 ) ^ (N U N 2 )}. The 

idea is to approximate the equilibrium distribution of customers at 

node 1 in Q with the equilibrium distribution of customers in Qi eq . 

At first glance, the approximation technique described above may 

seem rather ad hoc. However, we can draw upon the development in 
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Section II to understand the inner workings of the test-bed model, 
and to show that, at least in principle, the approximation scheme is 
not altogether unreasonable. 

It is clear that we can think of node 1, in the test-bed network Q, as 
comprising two FIFO nodes with service rates that depend just on the 
joint state of these two nodes. Consider the subnetwork Qi of Q 
consisting only of node 1. Theorem 1 can now be invoked to determine 
the exact equivalent network Qi eq . Let 

(for every (n u n 2 ) ^ {N U N 2 )) (for every r E \1, 2\), 

p^i,n 2 ) _ ^ja customer of type r is in service at node 

2/(n Y ,n 2 ) customers in Q x j. 

Q ]eq is then an open network consisting of node 1. When there are rii 
customers of type 1 and n 2 customers of type 2 in Qi eq , then customers 
of type r (G {1, 2j) enter the network at the rate xj." 1,n2) where 

(for every ( ni ,n 2 ) *S (NuN,))^^ = p^'V 

When a customer finishes service in Q leq it leaves the network (Fig. 
5). 
The equilibrium distribution of customers in Q leq is exactly the same 
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as the equilibrium distribution of customers in Qi. Observe, though, 
that the form of Qi eq is the same as that of (J leq , the difference lying 
in the state-dependent input rates. It is in this sense that the approx- 
imation scheme is reasonable. The idea now is to compare the exact 
state-dependent input rates, p&'^vi, with the approximate state- 
dependent rates, 

tN-fn,,^) m ( N r~ n r 



i.e., to compare pjg l,n * > with 



N r - n r 



Ni - Mi + N 2 - rh 
for all {n u ri2) % (N lt N 2 ) and for r E {1, 2}. 

3.3 Qualitative evaluation of the approximation 

In this section, we present a qualitative evaluation of the approxi- 
mation scheme as applied to the test-bed model. 

Observe that if the service rates for the two FIFO queues comprising 
node 1, in Q x , were not state dependent (in the priority scheme they 
are state dependent), then Qi would, in fact, be a product-form net- 
work. Theorem 2 would then lead us to conclude that Qieq and Q ieq 
were the same. Consider what happens if, in Q u "n is allowed to go to 
infinity. Then, effectively, the high-priority customers do not interfere 
with the low-priority customers at node 1. With u n - °°, the network 
becomes the one shown in Fig. 6, which is a product-form network. 
Thus according to our observation above, for values of v%\ that are 
large, compared to v x2 and v 2 , the approximation can be expected to 
yield very good results. 

In order to discover the situations in which the approximation can 
be expected to behave poorly, one needs to understand what aspects 
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of the exact network the approximation fails to capture. If Qi were a 
product-form network, then, given that (kuk 2 )(^(N\,N 2 )) customers 
were in node 2, all arrangements of customers within the node would 
be equally likely. As it stands, however, at node 1, priority 1 customers 
can preempt customers of priority 2. This suggests that given 
(ki,k 2 )(^Ni,N 2 ) customers in node 2, some arrangements of customers 
would be more likely than others. In fact, we conjecture that priority 
1 customers are more likely to be ahead of priority 2 customers, leading 
to the (conjectured) conclusion that 



(for every (0, 0) < (kM ^ WuNMdt*** > 



fe, 



k\ + k 2 



and 



p N-(*,.* 2 ) ^ 



k\ + k 2 



Thus $i eq uses smaller (resp. larger) state-dependent input rates for 
type 1 (resp. type 2) customers than the exact equivalent Qieq. This 
idea is suggestive, but it is difficult to draw any immediate conclusions 
from this conjecture as to the relationship between exact and approx- 
imate performance measures of the network. 

Another approach to discovering the direction in which the approx- 
imation can be expected to err is to observe that if node 2 in Q is 
replaced by a processor-sharing node, with class-independent service 
rate v 2 , then Qi^ becomes the exact equivalent of Qi (cf. Fig. 7). (This 
follows because when node 2 is processor sharing, if (k u k 2 )(^ (N U N 2 )) 
customers are present at node 2, then the rate of flow of class r(G 
jl, 2j) customers into node 1 is (k r /(ki + k 2 ))v 2 .) To fix ideas consider 
the case N x = n(^l) and N 2 = 1. The throughput of the class 2 
customer is simply the reciprocal of the mean successive passage times 
of the (single) class 2 customer through the point X (cf. Figs. 2 and 



(/V,. A/ 2 ) CUSTOMERS 




Fig. 7 — Q with node 2 processor sharing. 
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7). In either network, when the class 2 customer crosses the point X 
to enter node 2, it finds all the class 1 customers receiving service at 
this node. In the original network, since node 2 is FIFO, the class 2 
customer will have to wait for full service completion of the n class 1 
customers before it can leave node 2 (and subsequently, at some future 
time instant, cycle back through X). Thus, the mean sojourn time of 
the type 2 customer, in node 2 of the original network, is (n + l)/v 2 . 
However, if node 2 is processor sharing, then on entering node 2, the 
customer of type 2 starts receiving service immediately at the rate v 2 / 
(n + 1), and continues to receive service at a rate v 2 /(k + 1) (0 =£ k «£ 
n) until it finally leaves. Thus, in this case, the sojourn time of the 
class 2 customer at node 2 is stochastically dominated by an exponen- 
tially distributed random variable with mean n + l/v 2 , and hence has 
a mean smaller than (n + l)/v 2 . We further expect, intuitively, that, 
after completing service at node 2, when the type 2 customer returns 
to node 1, it expects to find more type 1 customers at node 1 when 
node 2 is FIFO than when it is processor sharing. Given that the type 
2 customer finds k(0 «S k «£ n) type 1 customers on its arrival at node 
1, its sojourn time at node 1 does not depend on whether node 2 is 
FIFO or processor sharing, and increases with increasing k. Thus, we 
expect that the mean sojourn time of the type 2 customer at node 1 
will be larger if node 2 is FIFO than when it is processor sharing. The 
conclusion is that the mean passage time of the type 2 customer, 
through the point X, is larger in the original network than in the 
approximating network. 

To see the magnitude of the error this effect could cause, let iVi = 1 
and allow v x2 —>oo f v 2 —> oo, and v 2 /u i2 -* 0. Under these assumptions, 
in the original network, the class 2 customer will be blocked once (and 
only once) at node 1 each time it cycles through the point X. The 
average time it spends in the blocked condition is l/v u - The rest of 
the time in each cycle tends to 0. Hence, the mean response time for 
the class 2 customer in the original network is l/v u . If node 2 is 
replaced by a processor sharing node, then, each time the class 2 
customer cycles through X, it is blocked once (and only once) with 
probability 1/2. Hence, the mean response time for the class 2 customer 
in the approximating network is l/2i> n , thus yielding an error of 100 
percent. 

We do not yet have a simple but rigorous argument that would allow 
us to say conclusively that the approximation yields higher through- 
puts for low-priority customers. However, the arguments presented 
above do make the conclusion plausible. 

3.4 Numerical examples 

To examine how the approximation works with specific examples, 
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we wrote a FORTRAN program to solve the equilibrium equations for 
Qieq using a simple recursive technique. 8 The program was somewhat 
more general, in that it could accept arbitrary state-dependent input 
rates and output rates. Thus, the same program could be used to solve 
the network exactly, if it were given the exact values of p[ 2 i,rt2) and 
p£»i."s) f or the various feasible (n u n 2 ). 

It is not hard to calculate exactly the probabilities p 2 V ,r%2) and pj§ 1,na) 
for some simple cases. Of course for (ni,n 2 ) # (Ni,N 2 ), p 22 1,n2) = 1 — 
p£i 1,n2) . Consider, for the purpose of illustration, the case Ni = 1, N 2 = 
1. Note that the state of the network Q is completely described by the 
state of node 2. The epochs of entry into the state S 2 = (12) are 
renewal epochs. The next state is, inevitably, S 2 = (2). The next state 
is S 2 = (21) with probability vu/(v\\ + "2) and S 2 = (0) with probability 
vz/(v\\ + v 2 ). Because of the preemptive discipline, the next state to be 
entered in the set ((12), (21)) will be S 2 = (12), thus completing a 
renewal cycle. Since the expected holding time in each state in {(12), 
(21) j is l/v 2 , therefore 

1 



.(0,0) _ 

P21 



v-> 



— + Vn i. 1 4. " n 



V 2 V\\ + V 2 V 2 V\\ + v 2 

and, of course, 

p 2 V 0) = and p 2 r = 1. 

In Table I we list the exact expressions for p 2 V'" 2 \ for all {n u n 2 ) *S 
(N U N 2 ), for some values of (N U N 2 ). These were computed in the same 
fashion as in the above example. 

In Tables 11(a), (b), and (c), we give several numerical examples of 
exact and approximate solutions of the test-bed network. The exact 
solutions and the approximate solutions were obtained using the 
FORTRAN program described above. The program yields the equi- 
librium joint-probability distribution of queue lengths at node 1. In 
Tables 11(a), (b), and (c), we display these joint probabilities and the 
node 1 utilizations. 

The following observations are immediate and summarize our con- 
clusions regarding the performance of the approximation scheme when 
applied to the test-bed network. 

1. The numerical computations support our earlier observations 
that if vu is large, then the approximation can be expected to yield 
excellent results [cf. case (1) in each of Tables 11(a), (b), and (c)]. 

2. The low-priority utilizations are consistently higher, again sup- 
porting our earlier observations regarding the direction in which the 
approximation can be expected to err. 
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Table I— Exact expressions for p^ 1 '" 2 ' in the test-bed network Q. 
(cf. Thm. 1 and Fig. 2) 



(0,0) 



Comparison 
Ni N 2 (n u n*) Pn^' with p 21 ' ' 

1 1 (1,0) 

(0,1) 1 



1+ * 



I'll + I»2 



(2,0) 

(0,1), (1,1) 1 

(i,o) *i 

I'll "2 "11 



1 + 



I'll + "2 "11 + "2 "11 + "2 

1+ "" 



(0,0) ^^ ; M 

"11 . / "li 



"11 + "2 \ "11 + "2 

"2 / "11 



"11 + "2 \ "11 + "2 

2 (1,0), (1,1) 

(0,2) 1 

(o.i) 1 ai 

"11 "12 



1 + 



"12 + "2 "11 + "2 "12 + "2 

[ ' "ii "2 "a " I 

[l»ll + "2 "11 + "2 "11 + "2j 



(0,0) - £$ 

"11 + "2 "11 + "2 

its 



L "12 "2 + "llj 



I»U + "2 

*(P22 ,,n,) = 1 - Pir ,) if (»!.««) * (MJV,); p^ 1 = pIb^ = 0). 



3. When v n , v i2 and v 2 are comparable, then the approximation 
yields good results with errors in the utilizations in the neighborhood 
of 10 percent. 

4. Considerable errors in the low-priority utilizations can arise, 
however. Witness case 3 in each of the Tables 11(a), (b), and (c). With 
a very large low-priority service rate at node 1, the approximate low- 
priority utilization suffers from an error of 20 to 50 percent. 

5. For the range of examples studied, the equilibrium probabilities 
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Table II — Numerical comparisons of exact and approximate 
solutions of the test-bed network 









Equilibrium State Probabili- 
ties at Node 1 


Node 1 Utilizations 








Approxi- 






Approxi- 










Exact 


mate* 




Exact 


mate* 


Case 






State 


Proba- 


Proba- 




Utiliza- 


Utiliza- 


No. 


"ii 


"12 


"2 (ni,n 2 ) 


bility 


bility 


Class 


tion 


tion 








(a) 


(N, = l 


, N 2 = 1) 








1 


10 


1 


1 (0,0) 


0.614 


0.606 














(1,0) 


0.029 


0.028 


1 


0.064 


0.063 








(0,1) 


0.322 


0.331 


2 


0.322 


0.331 








(1,1) 


0.035 


0.036 








2 


1 


.5 


1 (0,0) 


0.231 


0.222 














(1,0) 


0.077 


0.056 


1 


0.462 


0.444 








(0,1) 


0.308 


0.333 


2 


0.308 


0.333 








(1,1) 


0.385 


0.389 








3 


1 


100 


2 (0,0) 


0.393 


0.488 














(1,0) 


0.196 


0.163 


1 


0.601 


0.504 








(0,1) 


0.0059 


0.0081 


2 


0.0059 


0.0081 








(1,1) 


0.405 


0.341 














(b) 


(N, = 2 


, N 2 = 1) 








l 


10 


1 


1 (0,0) 
(1,0) 


0.681 
0.044 


0.676 
0.043 














(2,0) 


0.002 


0.002 


1 


0.0754 


0.0748 








(0,1) 


0.244 


0.249 


2 


0.244 


0.249 








(1,1) 


0.027 


0.027 














(2,1) 


0.003 


0.003 








2 


1 


.5 


1 (0,0) 

(1,0) 


0.179 
0.083 


0.171 
0.065 














(2.0) 


0.024 


0.016 


1 


0.631 


0.618 








(0,1) 


0.191 


0.211 


2 


0.191 


0.211 








(1,1) 


0.250 


0.260 














(2,1) 


0.274 


0.276 








3 


1 


100 


2 (0,0) 
(1,0) 


0.170 
0.116 


0.217 
0.109 














(2,0) 


0.050 


0.036 


1 


0.827 


0.780 








(0,1) 


0.0022 


0.0033 


2 


0.0022 


0.0033 








(1,1) 


0.187 


0.188 














(2,1) 


0.474 


0.447 














(c) 


(N, = 1, 


N 2 = 2) 








1 


10 


1 


1 (0,0) 
(1,0) 


0.464 
0.015 


0.457 
0.014 














(0,1) 


0.316 


0.319 


1 


0.0495 


0.0487 








(1,1) 


0.016 


0.016 


2 


0.487 


0.494 








(0,2) 


0.170 


0.175 














(1,2) 


0.019 


0.019 








2 


1 


.5 


1 (0,0) 
(1,0) 


0.116 
0.028 


0.109 
0.108 














(0,1) 


0.177 


0.182 


1 


0.437 


0.418 








(1,1) 


0.070 


0.055 


2 


0.447 


0.473 








(0,2) 


0.270 


0.291 














(1.2) 


0.340 


0.346 








3 


1 


100 


2 (0,0) 
(1,0) 


0.473 

0.167 


0.584 
0.130 














(0,1) 


0.008 


0.010 


1 


0.517 


0.404 








(1,1) 


0.115 


0.090 


2 


0.0101 


0.0123 








(0,2) 


0.002 


0.002 














(1,2) 


0.235 


0.184 









* Based on Norton's Equivalent. 
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are never drastically wrong, and follow trends similar to the exact 
values. 

The test-bed model is hard to analyze exactly for population sizes 
larger than the ones considered. We have run detailed simulations of 
the test-bed model for larger population sizes, and the results of these 
simulations continue to support the qualitative observations we have 
made above. 

IV. CONCLUSIONS 

We have demonstrated an approach for systematically analyzing 
the equivalent flow approximation. Our investigations have (1) re- 
vealed the conceptual basis for the approximation scheme, and (2) led 
to an understanding of the reasons for, and directions of, the errors 
that such an approximation scheme could introduce when applied to 
a class of prioritized central-server models. The approximation as 
described in the paper is of more general applicability, and much work 
remains to be done to discover its validity (accuracy and computational 
tractability) for more complicated, nonproduct-form networks. Our 
work, we think, provides the theoretical understanding and motivation 
for pursuing more detailed investigations. 
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APPENDIX 

Proof of Theorem 2: Index the nodes of Q£ in the same order in which 
they were indexed in Q. We need the following lemma. 
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Lemma: Let T {r) denote the class r routing probability matrix for the 
network Q 2 . Partition P (r) as follows: 

Mi M a 



p(r) _ 



M, 



P\ r l 



p(r) 



PS. 



p(r) 



22 



M 2 | i 21 

77ien 
(l)rGR 2 =» T ( '> = Pg + P$?[/ - PjJrPg 
(fi) // A (r) sota A (r) P (r) = A (r) £/w?n, partitioning A (r) as 
[Ai rt A 2 ' r 

. M x M 2 J 
(a) ifr£R 2 then A 2 r) = 
(fc) i/r G R 2 t/ien A 2 r) = A 2 r) T (r) . 
Proo/ o/ Lemma: Conclusion 1 follows readily from the fact that, for 
each r G R, P (r) is the transition probability matrix of a finite Markov 
chain with a single, positive communication class that has a nonempty 
intersection with {Mi + 1, • • • , M). For details see Ref. 7. 

Conclusion 2 follows directly from Conclusion 1. □ 

Returning to the proof of Theorem 2, we let ir. F& — * (0, 1) be the 
equilibrium distribution of the state process of network Q* It is now 
well known (cf. Ref. 9) that ir( ■ ) is of the form 

i M 

tt(S) = - n fi(s% 

U ,=i 

where G is a normalization constant and, for each i € {1, • ■ • , M), fc 
depends only N{S l ), v t and (for every r, (1 =£ r ^ R))\\ r \ where A (r) = 
(X< r) , • • • , \m) is any solution of A (r) P <r) = A (r) . 

Hence, (for every Si G F&) (for every i, r, (Mi + 1 < i < M), r G 
R 2 ) 

II fj(S j ) 



P*- 



|S':6'eFR,(i". • • • 1 .S'*')=6- 1 ,S'(l)=r| j-1 
M 

i n fj(s j ) 

|6- : SeF8,(.S 1 ,--,S v ')=S 1 |7'=l 

„ s n fjis$ 

IS»:Sj6F?nS^. w2(Si ,^i( 1 )=r| /=M,+1 

jl n fAsi) 

|.s- J :.s'2GFyn.s^L v2(Sl) l /=Af,+l 



*For notation see Section 2.2. 
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(where F^ H SjJl_ N 2 (Si) is the set of feasible states of Q 2 when 

N 2 - N 2 (Si) customers are in Qi), 

which, using the above lemma and the fact that the equilibrium 
distribution of the Q' 2 network state process is still product form, 

_ t N 2 -N 2 (S,) 
Kir 

Remarks: Some care is needed in asserting the last equality in the 
case where there are classes r t and r 2 , such that the submatrices of 
the communicating classes under T (ri) and T (r2) are the same permu- 
tation matrices (i.e., members of classes r x and r 2 cannot overtake 
each other). In this case the equality follows because for each N', flj 
and £*' are independent of the order in which members of these classes 
circulate in the network Q 2 . ED 
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